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Dynamics of localized states in extended supersymmetric quantum mechanics with 

multi-well potentials 

V.P. Berezovojj and M.I. Konchatnij 
A.I. Akhtezer Institute of Theoretical Physics, NSC "KIPT", 61108 Kharkov, Ukraine 

In this paper we propose a self-consistent approach to the description of temporal dynamics 
of localized states. This approach is based on exactly solvable quantum mechanical models with 
multi-well potentials and their propagators. States of Hamiltonians with multi-well potentials form 
a suitable basis for the expansion of wave packets with different shapes and localization degrees. We 
also consider properties of the tunneling wave packets, taking into account all states of Hamiltonians 
with symmetric and asymmetric potentials, as well as their dependence on the degree of localization 
and deformations of potentials. The study of the dynamics of initially localized states shows that 
application of the two-state approximation for the description of tunneling is considerably limited, 
especially for systems, which have several states in the under-barrier region, as for example in 
modern superconducting quantum interference devices and traps for cold atoms. 



I. INTRODUCTION 



Being one of the most exciting manifestations of the wave properties of particles, tunneling seemed to be a para- 
doxical prediction of quantum mechanics. Historically, the theory of the nuclear a-decay was the first study of tunnel 
"^^ ■ transitions pj], caused by bonding between nuclear resonances and continuous spectrum. Explanation of the inter- 
' \ molecular rearrangement of the ammonium spectrum [2| due to the tunnel splitting of vibrational spectra, opened 
iT^. the era of study of tunneling in multi-well potentials (3J. This problem became even more important in view of 
ri interpretation of the observed tunneling phenomena in condensed matter [4, 5] , and in particular in the Josephson 

junctions [6|. Currently, investigations of tunneling processes are mostly related to the study of the Bose-Einstein 
condensate in different types of traps 0, Q • 
►^ , Theoretical analysis of tunneling in double- well potentials is mostly performed within the two-mode approximation 

fy^ ■ Q. Important characteristics of this approach are the difference of energies of the ground state and of the first 
p^ \ excited state {/S. — Ei — Eq), as well as their wave functions. Values of A define, in particular, revival times of the 
l/~j ■ wave packets. This approximation explains general properties of tunneling, but it is unable to give interpretations 
CN ' of many subtle effects. Analysis of processes in multi-well potentials is complicated, since most of the models use 
r — , phenomenological or piecewise potentials (e.g., building from the rectangular wells and barriers, or from parabola), 
f^ '■ which are far from the real potentials. Spectra and wave functions in such potentials are unknown, that supposes the 
numerical analysis of their properties in consequent studies. It is important to note that there exist exactly solvable 
models with multi-well potentials |9l4l2l|. which could be used for the description of tunneling processes. 

Going beyond the two-mode approximation reveals interesting features of the dynamics of localized states. In 
particular, it leads to the substantial coherence breaking, as in symmetric [l3| as well as in asymmetric double- 
well potentials [IJI - llq . Moreover, characteristics of tunneling strongly depend on the shape of the potential, and 
^ . demonstrate a non-regular behavior which realized in bouncing increasing the probability of the localization of the 
wave packet (WP) in certain well. These features become even be more striking upon the initial localized state 
squeezing. 

In QM the time evolution of wave packets is described by propagators [2(j\, in which the contribution of the 
entire spectrum of the considered Hamiltonian is taken into account. However, only few models, mostly with one- 
well potentials, are known to have analytical expressions for propagators. To build propagators for exactly solvable 
quantum mechanical models with double-well potentials one could use the same approach as in |10l . llTl |l8| and 
construct new propagators from the known ones. 

The aim of the present paper is to describe the dynamics of localized states in multi-well potentials within the self- 
consistent approach. We use multi-well potentials obtained in the framework of iV = 4 SUSY QM UM and describe 
the dynamics of wave packets with corresponding propagators, calculated by use of the approach of |10l. llTI, Il8l|. 

The structure of this paper is as follows. In Section 1 we briefly discuss the construction of exactly solvable models in 
the framework of A'^ = 4 SUSY QM [12[ and give expressions for Hamiltonians with both symmetric and asymmetric 
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potentials, and wave functions, obtained from the initial model of harmonic oscillator (HO). In Section 2 we obtain 
the expressions for propagators in these models using the approach of [lO, [l3, 13 • Using the Hamiltonian of HO as 
the initial one, we obtain the explicit expressions for propagators in models with multi-well potentials. In Section 3 
we use the obtained potentials and propagators for the non-perturbative analysis of the dynamics of localized states 
in both symmetric and asymmetric potentials, and for different types of wave packets. In Section 4 we discuss some 
problems related to the considered problems and future developments. 

II. iV = 4 SUSY QM AND MULTI-WELL POTENTIALS 

To construct multi-well potentials in A^ = 4 SUSY QM [l^l we add states with energy e below the ground state 
energy Eq of the initial Hamiltonian Hq (we assume that Hq possesses only the discrete spectrum). Multi-well structure 
of the obtained potentials becomes more striking when -^^-g — - <C 1. The super-Hamiltonian of iV = 4 SUSY QM 

consists of three non-trivial Hamiltonians ^2j| H^ = H^ = Hq—e, HZ and H^. Spectra of the latter two Hamiltonians 
have the additional state below the ground state of the initial Hamiltonian, when iV = 4 SUSY is exact, while other 
states coincide with those of Hq — s. HZ and its wave functions are related to H^ = Hq — e and its initial wave 
functions in the following way: 

_ _ (f 

H_ = H^ - -^\n{ipi{x,e) + ip2ix,e)), 

w (x,E) — — = , 

v/2(i?. - s) ^{x,s,l) 

i:^{x,E = 0) = — — , (p{x,e,l) = <^i(x,e) +<f2{x,e) 

ip{x,e,l) 

-2 2W{ipi,ip2} ., . ipi{x,e) -^2{x,e) 



^-^- ..,:':r^.;^. ,^ . A(x,s 



A(+oo,£) - A(-cx),£)' ' Lpi{x,e) + Lp2{x,e) 

Here </?i(a;,£), i — 1,2 are two linear independent solutions to the auxiliary equation Hoip{x) = £1^9(2;). They are 
non-negative and have the following asymptotic behavior: ipi{x) — >■ -t-00 {tp2ix) — >■ 0) under x -> —00, and ipi{x) — > 
{if2{x) — >■ +00)) under x — > +00. Here 'ip^{x, Ei) are normalized wave functions of ifoj W^i'/'i, ^2} is the Wronskian. 
Using the form-invariance of H^ and HZ 12] one could obtain similar expressions for i?^, ip'^: 



hX = H^ - —■ ln((^i(a;,£) -f A(£, A) (p2{x,e)) 



^+ix,E = 0) = ^ 



{ipi{x,e)+A{e,\)ip2ix,e)) 
1 W{i;+{x,Ei),cp,{x,e)+A{e,X)Mx,e)} ^^) 

^^^''' '^ ^2{E, - £) {Mx,e)+A{e,X)Mx,e)) 

^ A(oo,£)-A-(A + l)A(-oo,£) 
^'' A(oo,£) + A-(A + l) A(-oo,£)' 

where the parameter A is restricted to be A > —1, and the normalization constant is N^'^ = (1 + A) N'^. 

In what follows, we will use the Hamiltonian of HO as the initial Hamiltonian to consider the tunneling of wave 
packets. Thus, the solutions to the auxiliary equation are parabolic cylinder functions ipi{S,, e) = £'^(v2C)i ¥'2{£,, e) = 
D^{—^/2£^) and the Wronskian becomes W{ipi,'fi2} = rf-iA O' '^■'^sre r(— i^) is the gamma-function, ^ = ^/uTx, 

^{t £- 1) - D,[V20 + D,{-V20. (3) 

Note that A(±oo,£, A), entering A(£,A), are determined by the asymptotes of the solutions (/9i(^,£), iy92(Cj^)- The 
ground state wave function of HZ is 

7V-1 N-^ , 2J^ 

ip-(x,E = 0) = ^ ^— = — -, N-^ = -^ — -, 4) 

and wave functions of excited states are determined by ([T]), where ^Z^{x,En) — (^) * -r=^Dn{"/2£^), n ~ 0,1,2... 
are the wave functions of HO. The potential of the Hamiltonian HZip,x) = uj (^5^ + (2 — ^) ~ -^'^'^fi^.T^i 1)) 




FIG. 1: Potentials and wave functions of Hamiltonians H_ and H'^. Left: (a) - potential U_ (^) [u] = 1, v = —0.02), wave 
functions of (b) - ground state, (c) - first excited and (d) - second excited states (d). Right: (a) - potential f/t(^,A + 1) 
{uj — 1, v — —0.02, A — —0.95) and corresponding wave functions. 

(see Fig. [T]) is defined by the symmetric com.bination of solutions ipi{x,e) + ip2{x,e), while the potential of H^ is 
defined by the asymmetric combination ipi{x,e) + (A + 1) ip2{x,e), and corresponds to the family of Hamiltonians 
with different values of to. Hence, the ground state wave function of Ht has the form of 



V'+(x,i? = 0) = 



N>:' 



n: 



{D,{V20 + (A + 1)D,{-V20) ^(C, £, A + 1) ' 



ivr^ = 



2{X + l)y/7ruJ 



(5) 



It should be noted that in terms of a dimensionless variable ^ the only way to vary the form of the potential is to 
vary < e < 1/2 and —1 < A. In the case of natural units x, additionally, the form of the potential (in particular, 
the position of local minima) can be changed by variation of w. 



III. TIME EVOLUTION OF STATES IN N = A SUSY QM 



In many cases tunneling is considered within the two-mode approximation. It allows to describe density oscillations 
and revival times for wave packets, but this approach fails, for example, in reproducing the dependence of the wave 
packet dynamics on the shape of the potential and the coherence breaking [IJ]. For instance, to correctly describe the 
dynamics of squeezed wave packets (WP) initially localized in one of the minima, one has to take into account not 
only the ground and first excited states, but also higher excited states. Contribution of these states increases with the 
squeezing of WP, and is also significant in potentials with low barriers, as for example in SQUIDs (superconducting 
quantum interference devices) and in cold atoms traps. 

/ 2R \ 1/4 / 

The time evolution of the Gaussian wave packet $(x) = I '^'^ J exp \—\{x 
parameter) initially localized in x = xq is determined by [201 



2^0) e^^] (here i? is the squeezing 



/+OC 00 
K{x,t;xo,{)) $(a;o) dxo, K {x, t; xq, 0) = y^ j'njx) V'n(xo) e" 
-00 , n 



iB„t 



(6) 



Here K {x,t;xQ,0) is the propagator, which sufficiently describes the dynamics of localized states in potentials of 
arbitrary complexity. 

Currently only a few exactly-solvable models with analytic expressions of the propagator K (x, t; xq, 0) are known 
[21| . Usually, their Hamiltonians have one- well potentials. The construction of new models and their propagators on 
the basis of propagators of exactly- solvable models with one-well potentials is proposed in [T^, [l3j 13 • Main ideas of 
this approach are outlined in ^Q, [13, [3 , so we will briefly discuss the procedure of the construction of propagators 
in. N = A SUSY QM, starting from exactly solvable model with the confinement potential. Let's denote propagators, 
corresponding to Hamiltonians H^^ oi N — A SUSY QM, as K^^ {x, t; xq, 0) (cr^ ~ ±). Using the form-invariance of 
A^ = 4 SUSY QM potentials, established in [l^l, the expression for the propagator K^ (x, t; xq, 0) could be obtained: 
K^ [x, t; Xq, 0) is related to the initial K^ {x, t; Xg, 0) of the exactly-solvable model as follows 



K 



+ {x,t;y,0) = -La:Ly 



+00 



dzK^ {x,t; z,0) Gj^(z,y,e) + 



N] 



-2 g-jrt 



<^(x,e,A) ip{y,e,k)' 



(7) 



Here L^ = [-j^ — ^ (^'^'ai ) ^^^ G^(z, y, e) is the Green function of the Schrodinger equation with energy e: 

2 
G^(x,y,e) ^ - „/r . r 1 {Mx,e) fr{y,e)0{y - x) + fi{y,e) fr{x,e)0{x - y)) . 

According to our notation, fi{x,e) = ip2{x,e), fr{x,e) = (pi{x,£). Acting with the operator Ly and simphfying the 
expression, the propagator becomes 



K+ix,t;y,0)^ T^TT ^- 



V 

A(£,A) / dz K^ {x,t;z,0) ip2{z,e)- 



dzKj^ (a;,i;z,0) ifi{z,e) 






^-2 g-^et (8) 



Lp{x,e,K)Lp{y,e,K)' 



Choosing the Hamihonian of HO as iJ^ , we get (^i(C, e) — D^{y/2£_), (p2{^, e) = D„{—\^^), A (e. A) = (1 + A), and 
K^ {x,t;y,0) becomes 21] : 

^+(-'^'^'")- ( 7.117 ) exp{^£^[(.^ + y^)cos.t-2.,]}, (9) 

where t = ^^^ + t, n € Nq,0 < t < —. The factor e~"*^2+") ensures correct behavior of the propagator for all 
values of time t and sewing of at ujt — Tin. Expressions dl]), ([S]), (O, ^ and © are the basic expressions to study 
the dynamics of localized states both in symmetric and asymmetric potentials. At the same time, this approach 
takes into account all the states, which form the localized state ^{x, 0). KZ {x, t; y, 0) could be obtained from (|S]) by 
substitution A (e. A) = 1 and corresponds to a symmetric potential. 

IV. DYNAMICS OF LOCALIZED STATES IN MULTI-WELL POTENTIALS 

In this paper we focus on the case, when only a few states of the Hamiltonian with the multi-well potential are 
located below the barrier. This is a common situation for different physical problems, both in atomic and in solid state 
physics. In general, the dynamics of the localized states can not be correctly described in the traditional framework 
of the tunnel splitting A = Ei — Eq, because the initially localized state can not be expanded as a superposition of 
the wave functions of the tunnel duplet. Higher excited states are essential in $(x, 0) and their contribution increases 
with increasing the localization of the wave packet. In contrast, the approach, considered in Section 3, allows one to 
study the dynamics of localized states taking into account all states of the exactly-solvable Hamiltonian, both with 
symmetric and asymmetric multi-well potentials. Since the spectra of HZ and H^ are identical, the tunnel splittings 
A = El — Eq are equal for symmetric and asymmetric potentials. 

We will study the dynamics of the localized states in potentials constructed from the potential of HO H])-®, with 
the propagator ([5]), ©• Their expressions in the dimensionless variables are 






Iri 

- ,1 ^ \ 1/2 



(A + 1) /' dCK-{^,T-C,0)D,{-^0 

•/ — oo 



(^(e,e,A+l)^(r;,e,A+l)' w' (10) 

X+(e,r;r^,0)=U-- exp--—[(C2 + 7?2) COST -2^7?] , 

\ 27r sm T j y2 sm 1 J 

where T = Tm + r, n Cz N, < r < tt. The propagator KZ {S,, T\ rj, 0) describing the dynamics of wave packets in the 
symmetric potential can be obtained from (llOp with A = 0. Thus, the expression ^ in the dimensionless variables 
reads 

/• + 00 / 2fl\l/4 2 

$(e, r) = y_ dCif+(C,r;C,0) $(C,0) , $(C, 0)=(^— j e-^, a^ =e-'''. (11) 

These relations allow one to obtain the form of the wave packet as a function of time and spatial coordinates. 



Small squeezed states. 

Let's note that study of the dynamics of the locahzed states (|lip can not be performed in the two-mode approxi- 
mation, even when R — Q. A satisfactory approximation of the initiahy locahzed state $(^, 0) is achieved with taking 
into account eight states of the Hamiltonian HZ (see line 1 in table 1). At the same time, the expansion of the 
initial state $ {£,, 0) over states of HZ (-ff+) is not required to compute |$ (^, T)\ when using the exact propagator. 

J2n=0 Cn V(Z\1{^, En) ex.p{-iEnT) 



"(-) + 



We will compare the results of the exact calculation using ([TT|) with |$(^, r)| 

to demonstrate the effectiveness of the basis of HZ {H^) in the considered problem. Figl^] shows |$(^,r)| for the 
same, as on Figdl potentials UZ{0 and U^{^,X), and for the value of squeezing parameter R — 0.35. The latter 
corresponds to the weak localization. Initially (at T = 0) the wave packet is localized in the left local minimum and 
has the energy (iJ$ — {^\HZ{p,x) |$)) E,s> = 0.16 for the symmetric potential, and i?$ = 0.185 for the asymmetric 
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-0.664 


0.018 


0.017 


0.048 


-0.146 


0.203 


-0.184 


0.110 


-0.031 


A = 0, iJ = 0.35 


0.682 


-0.692 


0.135 


-0.082 


0.056 


-0.069 


0.094 


-0.090 


0.044 


0.015 


A = -0.95, 1 


0.208 


-0.945 


0.179 


-0.105 


0.065 


-0.066 


0.077 


-0.060 


0.013 


0.029 


A = -0.95, r 


0.941 


0.213 


0.038 


0.011 


-0.007 


0.011 


0.066 


0.128 


0.153 


0.115 



TABLE I: Coefficients of the expansion of the initial wave packet (|lip . 




FIG. 2: Time dependence of wave packet [^(^jr)! with R = 0.35, initially localized in the left minimum of (a) 
UZ (C) {(.I = -2.29), (b) - potential [/+ (^, A) with A = -0.95 {^i = -2.153). 



potential 



Approximation of the initial WP by states of HZ and H^ is mainly determined by few lower states (see Table H]) 
which indicate the suitability of the basis. Nevertheless, despite a small contribution of higher states, they determine 
fine details of time evolution of WP, as we will show later. These details include specific beats caused by the 
interference of excited states. When the asymmetric potential is used, the expansion of the initial packet contains 
more terms, than in the symmetric potential (see Table |l|, and at the same time the contribution of states with ri > 9 
is comparable to the contribution of low-lying states. 

Thus, the temporal dynamics of WP (Fig. [2^) demonstrates slow tunnel transition of under-barrier states and 
fast oscillations of over-barrier states, which have higher intensity in compare to the symmetric case. In the case 
of a symmetric potential the evolution of I*!" {S,,T)\ (Fig. 2a) has the striking oscillatory nature: the portion of the 
initial $ (^, 0) which tunnels to the right minimum of UZ {S,) is quite large and reaches its maximum at T ~ — j^ 



\^ re 



300), while at T ~ Trev the WP is completely restored at the left minimum of U_ (^). At the 



(Ei-Ea) 

same time, the contribution of higher excited states to |$ {£,,T)\ is relatively small and leads only to small "beats". 
Completely different dynamics is observed, when the WP is initially localized in one of the minima of C/+ (C, A) (Fig. 
Wp)- The portion of the WP which tunnels to another local minimum is small, since in the left local minimum of 
C/t(^, A) (at ^/ = —2.153) the largest contribution to $ (^, 0) comes from the first excited state of Ht with the wave 
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300 




FIG. 3: Time dependence of |$(C,r)| for distributed initial packet ([12]) in potential U^ {(,, A) with A = -0.95. 

function being very small in the right well. It means, that this state contributes only a small portion to the tunnel 
transition amplitude. The similar situation is observed, when the wave packet is localized at the right local minimum 
{^r = 2.755) of f/+ (^, A). In this case the ground state wave function of H^, which is small in the left well, contributes 
mostly to $ (CjO)- ^^ some sense the wave packet is trapped within the initial well. The mechanism of such partial 
trapping of WP is simple: if one of the under-barrier states mostly contributes to $ (^,0), then its wave function is 
small in another well, that means that it has a small contribution to the tunnel transition to another well. Other 
under-barrier states give small contributions to tunneling due to their small portion in the formation of $(^,0). 
Nevertheless, the contribution of over-barrier states is larger than that in the symmetric potential, and it leads to 
large beats. 

The phenomenon of partial trapping is more obvious when the initial wave packet is uniformly distributed among 
the local minima (^^ , ^j?) of the asymmetric potential U^ {£_, A), e.g. 



H^,o) 



-2\ 1/4 

a \ e 



47r 









yo 



^2 = e-2fl. 
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(12) 



The dynamics of tunneling of the initial $ (^, 0) (|12p differs in different local minima of U^ {£_, A) (Fig. 3). 

Fig. [3] shows that the portion of WP ([T2| behaves differently in different wells of U^ (^, A). In the left minimum, 
for example, |$(^,T)| clearly oscillates, but at the same time the dynamics of |$(^,T)| in the right minimum is 
more complicated. To analyze it we give the time dependence of the square of the wave packet in minima of the left 
(6 = -2.153) and the right (^r = 2.755) wells (Fig. g]). 

It's important to note that in the left well wave packet oscillates and is completely restored after the time interval 
Trev — (E "^-E ) ~ '^^^ ' ^^ ^^® Same time, the fraction of the wave packet in the right well increases at < T < 150 
while |$(^,r)| is completely changed due to the tunneling from the left well. Moreover, when the value of time is 
close to T '-^ 150, the strong squeezing of the packet occurs, and under increasing of T |$(^, T)\ decreases and reaches 
its initial value. It looks like partial "confining" of the portion of WP inside the right well. Thus, in one of the 
minima (the left one) the tunneling dynamics possesses the oscillatory nature while in another minimum the partial 
"trapping" of a part of WP occurs. 

Another important characteristic of the WP tunneling dynamics is the probability to find wave packet in a certain 
well 



|$(2;,T)| dx. Fig. [5]shows P;(r)(2^) for the potential C/^(^, A) and the initial localized state 
(1121) . Pi(^r){T) is a quantitative characteristics of the localized state dynamics, and it can be revealed that a packet 



Plir){T)=J, 



xGl{r)well 




10 250 300 




FIG. 4: Time dependence of |$(C,T)| with R = 0.35 in (a) ~ left well (^; = -2.153), (b) - right weU (^r = 2.755). 
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FIG. 5: Pt(r){T) for the initial packet $(C,0) CU with R = 0.35. 

with the same distribution evolves in different local minima of U^{(,, A) in different way. The probability to find the 
packet in the right well initially increases due to the tunnel transitions from the left well, and then returns to the 
initial value during the time Trev — i^ '^-e ) ~ ^^^- Meanwhile, a portion of the WP in the right well, which tunnels 

to the left well C/^(^, A), is much smaller than in the reverse transitions, and the dynamics has typical oscillatory 
nature in the left well . This indicates that the partial trapping of the packet occurs in the deeper well, and typical 
oscillatory dynamics is kept in the left well t/+ (■Cj A). 

The initial distribution ([T2|) could be sufficiently approximated by ten states of the Hamiltonian iJ^ (See table HI]) , 
thus the proposed basis is effective enough for the description of localized states. Distributions like (fT2|) are typical in 
the study of the macroscopic tunneling of the Bose-Einstein condensate. Let's note that |<&(^, T) 

to ([TT1) coincides, with a good accuracy, with |$(^,r)| — ^^"q"' c„ ?/;S ^ (^,_E„) exp(— ii?„r) 
according to Table HIl 
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A = 


0.976 





0.165 








0.110 





0.058 





-0.061 


A = -0.95 


0.813 


-0.517 


0.154 


-0.067 


0.042 


-0.039 


0.101 


0.049 


0.117 


0.101 



TABLE II: Coefficients of the expansion of initial wave packet (|12ll . 



Large squeezed states. 

When the degree of the wave packet localization increases, the number of states, which significantly contribute to 
$ (f , 0), also increases. Let's consider the case when the initially localized state is located at one of the minima of the 
symmetric potential and in the deeper minimum of the asymmetric potential. When i? = 1.5 the WP with the center 
at .Jo = —2.29 is sufficiently approximated by twenty states of the Hamiltonian HZ (for a symmetric potential). The 
similar localized state with the center located at ^o = 2.755 needs twenty five states of Ht for accurate approximation 




FIG. 6: Time dependence of |$(^,T)| for initial packet with R = 1.5, localized in (a) - left minimum of potential U_ (^) 
(Ci = -2.29), (b) - in right minimum of [/+ (C, A), A = -0.95 {£,i = -2.153, £,,■ = 2.755). 
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FIG. 7: P;(,.)(r) for the initial wave packet "1> (^,0) with R = 1.5 initially localized in (a) - left minimum of the symmetric 
potential, (b) - in right minimum of the asymmetric potential. 



(for an asymmetric potential). Though the ground and the first excited states give the leading contribution to the 
expansion of wave packets, the contribution of higher states is still significant. It can be directly seen from the 
temporal dynamics of the initially localized states (Fig. [5]). 

The dynamics in the symmetric potential is complex, and the process of the barrier crossing can not be even called 
tunneling. Such a dynamics could be observed in the case of the coherent tunneling breaking in periodically driven 
systems [23,24], e.g. the contribution of the over-barrier states completely mimics the tunneling of the under-barrier 
states. Nevertheless, the revival time of wave packets is still Trev = i^ "^^^ \ ~ 300 according to the two-mode 
approximation. For the asymmetric potential the portion of the wave packet outside of the deeper minimum is small 
and its structure is complex enough. At the same time, the wave packet in the global minimum squeezes at time 
scales T ~ ^^ and restores at T = T^ev ■ 

This phenomena could be illustrated by the probability to find the particle in a certain well Pi(r) {T) — 
Ixei(r)weU l*^(^'^)l '^^ (Fig- El)- Large amplitudes of beats gives the evidence for the substantial contribution from 
the over-barrier states to the dynamics of wave packets. As it was mentioned above, during the evolution in the 
symmetric potential, a substantial part of the wave packet leaves the left well and return to the initial value at 
the time scale Trev Meanwhile, the portion of the wave packet in the right well of the asymmetric potential varies 
insignificantly, that tells on the partial trapping of the wave packet. 



V. CONCLUSIONS 

In the present paper we propose the approach to study of the dynamics of initiahy locahzed states, which is 
based on the exactly solvable quantum mechanical problems with multi-well potentials and on the corresponding 
exact propagators. Using the Hamiltonian of the harmonic oscillator as a basis one, we obtain, in frameworks of 
A^ = 4 SUSY QM, new Hamiltonians with multi-well potentials, both symmetric and asymmetric, together with 
the corresponding propagators. The study of the dynamics of the initially localized states demonstrates that the 
application of the two-mode approximation to the description of tunneling is very restricted, especially for systems 
with only few states in the under-barrier region. Such condition is typical for modern superconducting quantum 
interference devices and cold atoms traps. So, even the non-squeezed wave packet (i? = 0) can not be adequately 
approximated by wave functions of the ground and the first excited states. 

It is important to note that the states of Hamiltonians HZ and H^ are well-suited as the basic states to 
expand the localized states $(^,0). Usually, ten states are enough for a good approximation of <& (^, 0) in 
the sufhciently wide range of the squeezing parameter R. This can be confirmed by a good agreement between 

|$(^,T)| = ^^^=0" '^ni'(_' :){£,, En) exp{—iEnT) and the results of calculations by use of the exact propagator. In 

contrast to |13l4l6l |. where tunneling is compared in symmetric and asymmetric potentials, the spectra of Hamiltonians 
Hz and H^ are similar in our approach. Moreover, the shape of the potentials may be varied by the variation of e 
and A. 

The dynamics of WP contains slow tunnelin g of under-barrier states and fast beats, due to over-barrier states, and 
has a regular character, in contrast to [1^ [l^, [i6|. It happens since the over-barrier states have the spectrum of the 
initial Hamiltonian of HO (i.e. the equidistant spectrum) and the states get interfered with each other. Beats become 
smoother when a dissipation is taken into account [25| . The revival time of the wave packet is equal to T^ev = /^ ^^f n , 
that coincides with revival time predicted by the two-mode approximation. When the squeezing parameter is low 
[R = 0.35) the amplitude of beats is relatively small, since the number of excited states, contributing to WP is low. In 
the symmetric potential a part of the WP, which is formed by tunneling of the under-barrier states is smooth enough. 

In the asymmetric potential the dynamics of WP, initially localized in one of the minima of f/^T (^, A), has a number 
of distinctive features. In particular, the partial trapping of WP in the initial well and the suppression of tunneling to 
another well is observed. This phenomenon occurs independently on the well initially containing the WP. When the 
initial state is uniformly distributed between both wells of U^ (^, A), this phenomenon is observed for the deeper well, 
thus the tunneling rate from the deeper well to the other well is smaller than the rate of reverse transitions. When 
the squeezing parameter R increases, beats, caused by the over-barrier states, increase and mimic the contribution of 
slow tunneling of the under-barrier states. In some sense it could be considered as the destruction of tunneling of the 
initially highly localized wave packet. 
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of the manuscript and to A.J. Nurmagambetov for careful reading and correcting the manuscript. This research was 
supported in part by the Joint DFFD-RFBR Grant # F40.2/040. 
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